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IONOSPHERIC  MODELING 


E  and  F  Region  (100-500  km) 

The  previous  quarterly  report  reviewed  the  computer  code  that  was 
developed  to  generate  the  E  and  F  regions  of  the  polar  ionosphere. 

The  results  that  are  given  in  this  report  are  outputs  from  that  computer 
program. 

The  neutral  atmosohere  is  generated  from  the  atmospheric  model  of 
Jacchia  (1971).  Figure  1  shows  the  altitude  distribution  of  the  major 

neutral  constituents  based  on  an  exospheric  temperature  of  1000°K.  The 

*§ 

variations  of  the  atmospheric  constituents  with  respect  to  time  of  day, 
solar  cycle,  etc.,  can  be  fed  into  the  model  by  varying  the  exospheric 
temperature,  which  is  an  input  to  the  computer  program. 

There  are  two  processes  wnich  produce  Ionization  in  the  100-500  km 
height  region  of  the  polar  ionosphere;  solar  photon  absorption  and  ener¬ 
getic  particle  precipitation.  For  convenience  these  two  mechanisms  are 
further  subdivided  into  daytime  and  nighttime  photoionizction,  and  auroral 
arc  and  electron  "drizzle"  particle  precipitation. 

Figure  2  illustrates  the  photo  ion teat lor)  rates  versus  alt  tudes  for 
the  Jacchia  neutral  atmosphere.  The  daytime  production  rate  is  naturally 
a  function  of  latitude,  time  of  day  and  season;  the  conditiois  of  Figure  2a 
were  chosen  for  illustrative  purposes.  Figure  2b  records  the  ion  produc¬ 
tion  produced  during  the  nighttime  by  solar  radiation  scattered  to  the  dark 
side  of  the  ionosphere  by  the  earth's  hydrogen  and  helium  geocorona. 

This  proce; s  will  be  considered  constant  through  the  night  until  a  three- 
dimensional  photon  scattering  algorithm  Is  obtained  from  Robert  Meier 
of  the  Naval  Research  Laboratories  (1973). 
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Figure  2 


The  ionization  rate  distribution  produced  by  precipitating  particles 
is  given  in  Figure  3.  The  type  IBC-II  "standard"  aurora  (as  defined 
by  L.  R  Megill ,  Utah  State  University;  see  Ionospheric  Modeling  Newsletter 
for  April  10,  1973)  produces  ionization  near  100  km  at  an  order  of  magni¬ 
tude  higher  rate  than  the  noontime  solar  photon  flux.  The  "drizzle" 
electron  flux,  which  is  characteristic  of  the  polar  cap  *nd  perhaps  the 
whole  quiet  polar  ionosphere,  produces  the  ionization  shown  in  Figure  3b. 
Figure  4  is  a  plot  of  the  ion-pair  production  per  air  particle  per  incoming 
primary  electron  as  a  function  of  electron  energy.  This  figure  shows 
that  the  higher  energy  electrons  drive  deeper  Into  the  atmosphere  before 
reaching  their  peak  ionization  rate.  Since  the  atmosphere  has  an  expo¬ 
nential  density  variation  with  height,  the  total  production  rate  is  weighted 
toward  the  lower  altitudes  as  shown  in  Figure  3a  and  Figure  3b.  It  is 
obvious  that  the  time  and  spatial  distribution  of  precipitating  particles 
is  essential  in  the  determination  of  the  polar  ionospheric  distribution. 

A  copy  of  a  computer  code  which  would  include  the  transport  effects 
of  vertical  diffusion  was  found  to  have  been  superceded  by  a  new  code 
written  by  Mr.  Earl  Good  of  the  U.  S.  Air  Force  Cambridge  Research  Labora¬ 
tories  (1973).  This  new  code  has  been  promised  us  but  has  not  yet  arrived. 
This  will  put  our  diffusion  calculations  behind  schedule  but  we  still 
hope  to  have  a  diffusion  model  working  before  the  end  of  summer. 

D  Region  (30-1 00  km) 

Our  philosophy  for  modeling  the  very  complex  D  region  was  explained 
in  the  previous  quarterly  report.  A  subroutine  to  calculate  the  icn-electron 
pair  production  due  to  precipitating  protons  and  alpha  particles  has  bepn 
converted  f”om  the  (PAMATRX)  computed  program  discussed  by  Adams  and 
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Masley  (1965).  This  program  is  also  described  in  the  Ionospheric  Modeling 
Newsletter  of  May  10,  1973.  The  technique  for  calculating  the  D  region 
ionization  due  to  precipitating  electrons  is  the  same  as  that  *or  the 
E  and  F  regions  described  in  the  last  quarterly  report.  However,  the 
computer  subroutine  is  being  modified  slightly  to  be  compatible  with  the 
proton  subroutine  in  order  to  conserve  computer  time.  The  final  portion 
of  the  D  region  ionization  is  produced  by  solar  photons.  The  solar  flux 
has  two  major  effects  on  the  D  region  chemistry:  first,  the  ionization 
of  nitric  oxide  by  Lyman  o  and  second,  the  dissociation  of  heavy  molecules 
by  radiation  of  wavelengths  greater  than  1200  A.  We  have  found  that  if  the 
detailed  structure  of  the  02  absorption  cross  section  and  the  NO  ionization 
cross  sections  are  taken  into  account,  then  the  solar  flux  in  the  1030-1280  A 
band  may  contribute  as  much  N0+  ionization  as  the  Lyman  a  line  at  1215.7  A. 

This  will  be  a  very  significant  result  since  most  D  region  models  consider 
only  Lyman  a  production.  Spurred  on  by  these  results,  we  are  now  calcu¬ 
lating  the  ionization  rate  of  02(1Ag)  (an  excited  state  of  02)  which 
might  also  contribute  an  additional  25%  to  the  D  region  ionization  pro¬ 
duction.  Detailed  results  will  be  available  for  the  next  reprrt.  Modification 
of  the  E  and  F  region  photon  production  subroutines  is  being  made  to  accom¬ 
modate  these  latest  results  into  a  computer  code  for  the  D  region.  Once 
these  production  subroutines  are  complete,  in  the  next  quarter,  the  four-ion, 
steady-state  model  of  the  D  region  chemistry  can  be  run  to  produce  ionization 
profiles. 
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